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ABSTRACT 



We report the results of a Suzaku X-ray imaging study of NGC 625 1 , a nearby 
giant radio galaxy with intermediate FRI/II radio properties. Our pointing direction 
was centered on the 7-ray emission peak recently discovered with Fermi-LAT around 
the position of the north-west radio lobe 15 arcmin offset from the nucleus. After 
subtracting two "off-source" pointings adjacent to the radio lobe, and removing pos- 
sible contaminants in the XIS field of view, we found significant residual X-ray emis- 
sion most likely diffuse in nature. The spectrum of the excess X-ray emission is well 
fit by a power law with photon index T = 1.90 ±0.15 and a 0.5 — 8keV flux of 
4 X 10~^^ erg cm~^ s^^. We interpret this diffuse X-ray emission component as being 
due to inverse-Compton up-scattering of the cosmic microwave background photons 
by ultrarelativistic electrons within the lobe, with only a minor contribution from the 
beamed emission of the large-scale jet. Utilizing archival radio data for the source, 
we demonstrate by means of broad-band spectral modeling that the 7-ray flux of the 
Fermi-hAI source 2FGL J1629. 4-1-8236 may well be accounted for by the high-energy 
tail of the inverse-Compton continuum of the lobe. Thus, this claimed association of 
7-rays from the north-west lobe of NGC 6251, together with the recent Fermi-LAT 
imaging of the extended lobes of Centaurus A, indicates that particles may be effi- 
ciently (re-)accelerated up to ultrarelativistic energies within extended radio lobes of 
nearby radio galaxies in general. 
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8 Subject headings: acceleration of particles — radiation mechanisms: non-thermal — 

9 X-rays: galaxies — gamma rays: general — galaxies: jets — galaxies: individual 

10 (NGC6251) 



Introduction 



12 With the successful launch of the Fermi Gamma-ray Space Telescope, we now have an un- 

13 precedented opportunity to s tudy in detail the 7-rav emission from different types of astrophysical 

14 sources. The 2FGL catalog (|Abdo et alj|201 contains 1873 sources detected in the 100 MeV to 

15 100 GeV range by the Large Area Telescope (LAT) onboard the Fermi satellite during the first 24 

16 months of the science phase of the mission, which began on 2008 August 4. Among 1297 objects 

17 included in 2FGL which are firmly identified or reliably associated with counterparts of known 

18 source classes, more than half are blazars, i.e. active galactic nuclei (AGN) with radiative outputs 

19 dominated by a beamed jet emission. Non-blazar AGN constitute only a small fraction (~ 3%) of 

20 the identified 2FGL population, but nonetheless have already emerged as a new, important, and at 

21 the same time relati vely diverse class of 7-r ay emitters. The second catalog of AGN detected by 

22 Ferm/-LAT (2LAC: lAckermann et al.ll201 la|) contains ten radio galaxies, four narrow-line Seyfert 

23 1 galaxies, and several other non-blazar-typ e AGN with prominent starburst components (see also 



Abdo et al.ll2010d : lAckermann et al.ll201 Ibh 



25 Particularly noteworthy in this context is the Fermz-LAT detection of the giant lobes in the 

26 nearby radio galaxy Centaurus A, extending for about ~ 5° in the east-west direction, and ~ 9° in 

27 the north-south direction (~ 300 kpc x 600 kpc at the distance of 3.7Mpc: lAbdo et al.ll2010blJdh . 

28 The Fermz-LAT results regarding the Centaurus A lobes have profound implications for the general 

29 understanding of the evolution of radio galaxies, their interactions with the surrounding environ- 

30 ment, and the production of high-energy particles in the Universe. In particular, pro mpted by these 



31 recen t results from the Fermz'- LAT and also from the Pierre Auger Observatory ([Abraham et al 



32 



2OO7I : iMoskalenko et al.ll2009h . several authors have debated on the high-energy particle acceler- 



33 ation and produ ction of ultra high-energy cosmic rays within the giant lobes of C entaurus A and 

34 similar systems faardcastle et al.ll2009l : lo' Sullivan et aDl2009l : IPe'er & Loebll201ll) . 



35 



Besides Centaurus A, three other radio galaxies detected so far with Fermi-LNI could be 

36 promising targets for investigating the extended 7-ray emis sion components assoc i ated \y ith radio 

37 lobes. Prior to the launch of Fermi, ICheungI (|2007l) and iGeorganopoulos et al.l (|2008|) pointed 

38 out that the radio lobes of Fornax A might be detected at GeV energies and seen as extended 



See also the online version at http : / / f ermi .gsfc.nasa. gov/ssc/data/ access /I at/ 2yr_catalog/ 
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39 Structur es by the LAT bu t its associated 2FGL source is faint and currently appears to be a point- 

40 source (|Abdo et al. 2011 ). Similarly, the 7-ray emission associated with Centaurus B and reported 

41 recently in 2LAC may be dominated by the lobes, but due to a relatively small angular size of 

42 the source (~ 0.°2 in the north-south direction) careful analysis of the Fermi-hAT data is needed 

43 before drawing any robust conclusions (Katsuta et al. in prep). The giant radio ga laxy NGC 6251 

44 already appeared in the first-year Fermi-LAI catalog (IFGL; lAbdo et al.ll2010a ) and is another 

45 object worth attention due to the large angular extent of its radio structure (~ l.°2) which in 

46 principle could be resolved with Fermi-hAI above 1 GeV photon energies. Based on the one-year 

47 of accumulated LAT data, the NGC 6251 nucleus was included within the 95% LAT error ellipse 

48 of IFGL J163 5.4-I-8228, and its re ported 7-ray flux could be due to the unresolved nuclear jet of the 

49 radio galaxy (IMigliori et al.ll201 lb . Subsequently however, the position of the 2FGL counterpart 

50 of NGC 6251, namely 2FGL J1629.4-I-8236 (detected at 12 a level), is now shifted north-west with 

51 respect to the IFGL source position, toward the outer jet and radio lobe. We note that previously, 

52 the NGC 625 1 galaxy was proposed to be associated with the EGRET source 3EG J1621-I-8203 

53 dMukherjee et al.ll2002l) . 



54 NGC 62 5 1 is a nearby AGN hosting a s upermassive black hole with a mass, Mrh — (4 
8)x 



56 



0« Mq (iFerrarese & Fordl[l999l : lHoll2002 ). It is classified as an FR I radio galaxy (iLaing et alj 



9831) based on its jet morphology and monochromatic 178 MHz luminosity of 1.9 x 10 W Hz sr 



57 dWaggett et al.l ll977|) . alt hough the entire giant radi o structure resembles more FRII or interme- 

58 diate FRI/FRII systems dschoenmakers et allboOlk A t the distance of ~ 106 Mpc (z = 0.0247, 



59 conversion scale 0.491 kpc/arcsec; IWegner et al.l 



2003|) the linear extension of the lobes reads as 



60 ~ 2.1 Mpc, which is a larger physical size than th at covered by the o uter lobes of Centaurus A. 

61 The radio jet in NGC 6251, discussed in detail in (|Perley et al.lll984|) . is likely to be relativistic 



62 up to large distances from the core and observed at viewing angles 9j < 40 deg, as suggested by 

63 the jet-to-counterjet brightness asymmetry on both pc and kpc scales. It can be divided into (i) 

64 the bright inner region within the first ~ 120", (ii) a faint central part extending from ~ 120" 

65 up to ~ 180", (iii) the outer region flaring between ~ 180" and ~ 270", and (iv) a curved and 

66 low- surface brightn ess tail beyond ^ 270". The outer jet region ha s been detected at X-ray fre- 



67 quencies by ROSAT (IBirkinshaw & Worra' 



68 (|Sambruna et al.ll2004a : 



il993l:lMacketaL 



1997b ). XMM-Newton and Chandra 



Evans et al.ll2005|) . The origin of the detected jet X-ray emission is still 



69 an open question a nd different models have been discus sed for such, including thermal emission 

70 (IMack et al.lll997bb. beamed inyerse-Compton rad iation (|Sambruna et al.ll2004ab . and synchrotron 

71 process dSambruna et al. 2004a ; Evans et al. 2005). However, none of the previous X-ray observa- 

72 tions pointed to either the north-west or south-east lobe regions, separated by ~ 0°.3 and ~ 0°.7 

73 from the NGC 6251 nucleus, respectively. 

74 In this paper, we p resent the results of the newly conducted observations of NGC 625 1 with the 

75 Suzaku X-ray satellite ( Mitsuda et al.ll2007|) mainly in 201 1 April. NGC 6251 was first observed by 
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76 Suzaku in 2010 December (|Evans et al.ll201l|) . but these observations targeted strictly the nucleus 

77 of the radio galaxy. This time the pointing direction was intentionally shifted 15 arcmin in the 

78 north-west direction from the nucleus in order to study the 7-ray emission peak discovered with 

79 Fermi-LAT around the position of the NW radio lobe. In § 2, we describe the Suzaku observations 

80 and data reduction procedure. The results of the analysis are given in §3. The discussion and 

81 conclusions are presented in § 4 and § 5, respectively. 



2. Suzaku Observations and Data Reduction 



83 Suzaku is the ideal instrument for the intended study since it provides a low and stable No n 

84 X-ray Background (NXB), particularly for diffuse sources (IMitsuda et al.ll2007l : lTawa et al.ll2008|) . 

85 We have therefore observed the NW lobe in NGC 6251 for 34.2 ksec ("on" pointing). In addition, 

86 two short observations (18.8 ksec and 12.1 ksec) were made in the vicinity of the radio lobe in 

87 order to estimate the background flux ("off" pointings). Hereafter, we denote these three obser- 

88 vations as LOBE, BGDl, and BGD2, respectively, and the observation log is provided in Table 1. 

89 Both BGDl and BGD2 observations were made just after the LOBE observation so that the long- 

90 term background fluctuation as well as possible effects caused by the detector degradation can be 

91 minimized. In the case of BGD2, an additional short (11.1 ksec) exposure was obtained in 2011 

92 August, but these data were excluded from the analysis presented here to avoid possible artifacts 

93 related to the aforementioned problems. We have checked that the presented results do not change 

94 within statistical errors when the additional BGD2 data are included in the spectral analysis. 



95 The observations were made with X-ray Imaging Spectrometer (XIS; iKoyama et al.l |2007|) 

96 which consists of four CCD cameras each placed in the focal pla ne of the X-ray Telescope (XRT; 
Serlemitsos et aklboOTh . and with the Hard X-ray Detector (HXD: lKokubun et al.ll2007l : lTakahashi et al 



2007b . Each XIS covers a 18' x 18' region on the sky, which is sufficiently large to include the 



99 whole NW lobe within a single pointing imag e. On the other hand, s ince HXD is a non-imaging 



100 instrument with a large field of view 34' x 34' (ITakahashi et al 



101 between the hard X-ray emission from the nucleus (lEvans et al 



2007) and we cannot discriminate 



20111) and possible emission from 



102 the extended NW radio lobe in the HXD data, we only use the imaging and spectral results obtained 

103 with XIS in this paper. One of the XIS sensors (XIS 1) has a back-illuminated (BI) CCD, while the 

104 other three (XISs 0, 2, and 3) utilize front-illuminated (FI) CCDs. Because of an anomaly in 2006 

105 November, the operation of XIS 2 has been terminated and we use only three CCDs. The XIS was 

106 operated in the normal full-frame clocking mode with the 3 x 3 or 5 x 5 editing mode. 

107 In our Suzaku pointing, the relatively bright NGC 6251 nucleus was located well outside the 

108 field of view of the XIS. However, due to the r elatively large point sp read function of the XRT 



109 onboard Suzaku (half power diameter of ~ 2.0': lSerlemitsos et al.ll2007|) . the contamination of the 
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110 NGC 6251 nucleus may still, in principle, affect the exposure of the NW radio lobe: the pointing 

1 center of our LOBE observation is separated by only ~ 15' from the core. For this reason, we 

2 simulated the contamination effect using xissim. In the simulation, we assumed a point source 

3 at the position of the nucleus characterized by a power-law continuum with the photon index, T = 
4 1.87, and the normalization K = 6.79 x 10^^ photons keV^^ cm^^ s"^ at 1 keV (correspondin g 



5 to the 2 — lOkeV unabsorbed luminosity 2.8 x 10^^ erg s^^), as reported in (|Evans et al.ll201 II) . 

6 We found that the expected net count rate (the XIS count rate for the "on" pointing subtracted by 

7 the background count rate estimated from the "off" pointings) for the region used in the spectral 

8 analysis (see § 3.2 below) due to the core contamination is 4.11 x 10^^ s^^ only, whilst the total 

9 source net count rate is 6.40 x 10"^ s~^ Hence, we conclude that the leakage effect is at the 

120 level of about 0.64 %, and so that the stray light contamination from the NGC 6251 nucleus can be 

121 considered as negligible when analyzing the diffuse emission of the NW lobe. 

122 We analyzed the screened XIS data, reduced using the Suzaku software version 1.2. The 

123 screening was based on the following criteria: (1) only A5CA-grade 0, 2, 3, 4, 6 events were accu- 

124 mulate d, while hot and flickering pixels were removed from the XIS image using the sisclean 



125 script (iDay et al.l 119981) . (2) the time interval after the passage of South Atlantic Anomaly was 

126 greater than 60s, (3) the object was at least 5° and 20° above the rim of the Earth (ELV) during 

127 night and day, respectively. In addition, we also selected the data with a cutoff rigidity (COR) 

128 larger than 6 GV. In the reduction and the analysis of the Suzaku XIS data, HEADAS software 

129 version 6.11 and a calibration databases (CALDB) released on 2009 September 25 were used. 

130 The XIS cleaned event data-set was obtained in the combined 3x3 and 5x5 edit modes using 

131 xselect . We note that no spectral variability were observed in the XIS data and the imaging and 

132 spectral analyses are discussed in the next section. 



3. X-ray Analysis 



134 



3.1. Image Analysis 



135 



We extracted the XIS images within the photon energy range of 0.4 — lOkeV from only the 

136 two FI CCDs (XIS 0, XIS 3) because the BI CCD (XISl) has lower imaging quality. In the image 

137 analysis, we excluded calibration sources at the corner of th e CCD chips. The images of the NXB 



138 were obtained from the night Earth data using xi snxbgen (|Tawa et al. I boosl) . Since the exposure 

139 times for the original data were different from that of NXB, w e calculated the appro priate exposure- 

140 corrected original and NXB maps using xisexpmapgen dlshisaki et al.ll2007h . The corrected 

141 NXB images were next subtract ed from the correct ed original images. In addition, we simulated 

142 flat sky images using xissim (|lshisaki et al.ll2007h . and applied a vignetting correction. All the 
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143 images obtained with XISO and XIS3 were combined and rebinned by a factor of 4 (CCD pixel 



44 size 2 4 fim x24 fxm, so that 1024 x 1024 pixels cover an 18' x 18' region on the sky: lKoyama et al 



2007 '). Throughout these processes, we also performed vignetting correction for all the XIS images 



145L 

146 for LOBE, BGDl and BGD2. Finally, the images were smoothed with a Gaussian function with 

147 a = 0.28". 

148 The resulting X IS images are shown in Figure 1, with a 325 MHz radio contours of NGC 6251 



149 (iMack et al.lll997ah and the 95% LAT error circle of 2FGL J1629.4+8236 overlaid. When compar- 

150 ing the LOBE, BGDl, and BGD2 images, one can clearly see an excess X-ray emission possibly 

151 associated with the NW radio lobe. The apparent enhancements at the edges of the XIS images are 

152 artifacts due to an insufficient exposure at the CCD edges caused by small fluctuation of attitude 

153 during the observation. Interestingly, at the sensitivity limit of these Suzaku observations, no ex- 

154 cess X-ray emission was detected from the bright radio hotspot located at the edge of the lobe. In 

155 the context of the study by lHardcastle et al. I (|2004|) . since NGC 625 1 is a low power FR I/II source. 



156 the X-ray emission from the hotspot is expected to be dominated by the high energy energy tail of 

157 the synchrotron continuum rather than by the inverse-Compton component. Specifically, the radio 

158 flux of the hotspot of 0.1 Jy at 1.5 GHz allows us to estimate roughly the expected X-ray flux by 

159 extrapolating the observed hotspot spectrum from lower frequencies as ~ 2 nJy. A point source 

160 at this flux level is easily detectable by Chandra or XMM, but our Suzaku non-detection is unsur- 

161 prising due to the relatively poor angular resolution of the instrument combined with a relatively 



162 intense diffuse emission of the surrounding lobe. 



163 3.2. Spectral Analysis 

164 In the spectral analysis, we extracted photons from the arc-like source region in the LOBE 

165 pointing as shown in Figure 2. This 'arc-like region' we are referring to is a part of a concentric ring 

166 centered at the radio core position of NGC 625 1 . There we exclude d X-ray point sources detected 

167 by ROSAT (see the TABLEl described in iMukherjee et al.l |2002|) as background contaminants. 



168 assuming the source region radii 2' (a typical half -power diameter of the XRT, as described above). 

169 We note that the origin of these X-ray features is still under debate, and hence it remains possible 

170 that they represent hotspot and enhancement related to the NW lobe itself. In our analysis, however, 

171 we removed all such bright X-ray spots, since the main objective of the paper is to detect and to 

172 characterize the dijfuse X-ray emission component associated with the dijfuse radio structure of 

173 the lobe. Similarly, we defined arc-like background regions for both BGDl and BGD2 pointing 

174 using the same detector coordinates as chosen for the LOBE source region. We extracted the 

175 spectra from LOBE, BGDl and BGD2 source regions with xselect for each CCD (XIS 0, XIS 
176 1, XIS 3). Spectral analysis and model fitting were performed with xspec version 12.7.0. In order 
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177 to improve the statistics, the spectra of XISO and XIS3 were summed together using mathpha. 

178 Moreover, the spectra of BGDl and BGD2 were summed for the same reason in the same way. 

179 Finally, we made redistribution mat rix files (RMFs) an d ancillary response files (ARFs) using 

180 xisrmf gen and xissimarf gen dlshisaki et al.ll2007r) . respectively. 



181 Using these RMFs and ARFs, the corrected spectrum about energy response and the effec- 

182 tive area of XIS were obtained. Figure 3 shows the resulting (background subtracted) spectrum of 

183 the LOBE pointing within the energy range 0.4 — 7.5keV. The spectrum could be well fitted by 

184 a single power-law continuum with a photon index, T = 1.90 ± 0.15, moderated by the Galactic 

185 absorption only. The Gala ctic hydrogen column density was fixed as A^h = 5.54 x 10^° cm"^ 

186 (IDickey & Lockmanlll990|) . The value of x^/d-0.f= 20.54/14 indicates that this is a satisfac- 

187 tory model for the diffuse X-ray emission component of the NW lobe in NGC 6251. The details 

188 of the fitting results are summarized in Table 2. We note that, the spectrum is also equally well 

189 reproduced by a thin thermal plasma models (either bremsstrahlung and/or Raymond-Smith) mod- 

190 crated by Galactic absorption, with the temperature of kT > 3. 7 keV. This is rather high compared 



91 with those of thermal plasma of nearby elliptical galaxies (e.g. [Matsushita. Ohashi. & Makishima 



192 



2000) . Moreover, the implied thermal pressure, pth ~ 2 x 10 dyn cm ^, is more than four orders 



193 of magnitude larger than the minimum-energy non-thermal pressure of the lobes (see, section 4.2). 

194 We regard these parameters as unrealistic, and hence we conclude that the diffuse X-ray emission 



95 component detected at the position of the NW lobe is purely non-thermal in origin. 



196 



4. Discussion 



197 The analysis of Suzaku data for the NW lobe in NGC 6251 revealed the presence of X-ray 

198 emission, most likely of diffuse nature, well fitted by a power-law continuum with photon index 

199 r < 2 moderated by the Galactic absorption only. The 0.4 — 7.5 keV luminosity of the lobe corre- 

200 sponding to the Suzaku source extraction region reads as 5.4 x 10"^^ erg s~^. The detection of such a 

201 non-thermal X-ray emission component is interesting, but not exceptional, since several analogous 

202 dete ctions of the X-ray lobes i n either FRI or FR I I radio galaxies h ave been previously reported 

203 (see lKataoka & Stawaral2005l : ICroston et al.ll2005l : llsobe et al.ll201ll . and references therein). The 

204 relevance of our Suzaku observations of the NW lobe in NGC 6251 is due to the aforementioned 

205 possible (or even likely) association of the lobe with the Ferm/-LAT source 2FGL J1629.4-I-8326. 

206 Note that in the case of the only radio galaxy for which the lobe emission was robustly associ- 

207 ated with a 7-ray source so far, Centaurus A, no diffuse X-ray e mission component rel ated to the 

208 extended radio structure has been detected (see in this context. iHardcastle et al.ll2009h . In other 

209 words, our observations may potentially provide the very first case of the detected X-ray emission 

210 from the large-scale lobe bright in 7-rays. 
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211 Yet the above statement may be considered as premature, since the 95% LAT error circle 

212 of 2FGL J1629.4+8236 includes not only the extended lobe but also the bright 'outer jet' re- 

213 gion (see § 1 and Figure 1). This region, as mentio ned previously, has been resolved in X-rays 

214 before with ROS AT, Chandra, and XMM-Newton dMack et alJll997bl : Isambruna et all l2004al : 



Evans et al.ll2005[) . Regardless of the debated origin of the jet-related X-ray photons (synchrotron 



216 versus beamed inverse-Compton), large scale jets detected at keV photon energies are in general 

217 thought to be the sites of efficient acceleration of high-energy electrons, and as such, to be potential 

218 sources of 7-rays detec table by modern instruments like Fermi-LNI (see, e.g., Stawarz et al.ll2004 : 



219 



Sambruna et al. l2004bl) . Below we discuss in more detail the possible association of 2FGL J1629. 4-1-8236 



220 with the outer jet and with the NW lobe of NGC 6251, and demonstrate by means of broad-band 

221 spectral modeling that the 7-ray flux can best be accounted for by the high-energy tail of the 

222 inverse-Compton continuum of the lobe. 



4.1. Outer Jet 



224 For the spectral mod eling of the outer jet reg ion, we take the 1.4 GHz flux density, S'i.4ghz = 

225 0.5 ± 0.05 Jy measured by lSambruna et al.l (l2004a|) fr om a VLA map, an d assume the radio spectral 

226 index = 0.64 ± 0.05 within the 1—5 GHz range (iPerley et al.lll984|) . We also approximate the 

227 X-ray spectrum of this region as a power-law with the photon index F = 1.68 ± 0.1 3 within the 

228 0.4 — 7keV range and the monochromatic flux 5*1 kev = 4.7 ± 0.4 nJy, as given in Evans et al. 

229 J2005l) l^ We approximate the jet as a cylinder with a radius _Rj = 5 kpc and projected length 

230 £j = 40 kpc so that the emitting volume is Vj = TiRfij/ sin 9y 

231 We model the broad- band emission of the outer jet region in t he framework of the beamed 

232 inverse-Compton scenario (ITavecchio et al.ll2000l : ICelotti et al.ll200lb . with the radio emission be- 

233 ing due to the synchrotron radiation of non-thermal electrons distributed isotropically within a 

234 relativistic outflow, and the X-ray-to-7-ray continuum due to Comptonization of the cosmic mi- 

235 crowave background (CMB) photons by the same electron population ('beamed EC/CMB' model 

236 below). Anisotropy of the CMB radiation in the jet r est frame, as well as Klein-Nishina effects, 

237 are properly taken into account (iStawarz et al.ll2005|) . For simplicity, we assume homogeneous 

238 distributions of the radiating electrons and of the magnetic field within the whole outer jet region, 

239 which should be considered as a rough approximation only consi dering the observed asymmetric 

240 radio intensity and polarization transverse profiles of the outer jet (IPerley et al.lll984|) . 



^We note that slightly different X-ray spectral parameters for the outer jet region are claimed bv lSambruna et al. 
(l2004ah . namely Si kov — 2.3 ± 0.7nJy and F = 1.15 ± 0.3 8 within the 0.5 — 9 keV range. Our re-analysis of archival 
XMM-Newton data is however more consistent with that of lEvans et al.l (l2005b . 
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241 As shown in Figure4 (thin solid line), the 7-ray flux of 2FGL J1629.4+8236 is not well ac- 

242 counted for by the beamed EC/CMB model for the jet with the following 'best fit' parameters: the 

243 ratio between comoving electron and magnetic field energy densities r/cq = U'^/U'^ = 1.75, mag- 

244 netic field intensity B = 2.3 /iG, jet bulk Lorentz factor Fj = 3, jet viewing angle 9j = 17deg, and 

245 the electron energy distribution of a power-law form ^^(7) oc 7"'' x exp[— 7/7max] with s = 2.2, 
2t6 7min = 10 and 7max = 5 X 10^. With the above model parameters the observed ratio of the inverse- 

247 Compton and synchrotron peak luminosities reads as Ljc/Lgyn — (S/T^y x (U'^^^^/U^) ~ 25. Note 

248 however that the relatively small jet viewing angle implied by the model, 9^ = 17 deg, may be con- 

249 sidered as unlikely consideri ng the large projec ted size of the whole radio structure of NGC 625 1 

250 (see the related discussion in Evans et al.ll2005|) . On the other hand, the implied jet kinetic lumi- 

251 nosities stored in electrons, magnetic field, and cold protons (assuming equal comoving number 

252 densities of protons and electrons), namely Lc — 0.7 x lO''^ erg s^^. Lb — 0.4 x 10"''* erg s^^, and 
1.9 X 10^^ e rg s~^, giving the to tal jet power Lj ~ 2 x lO^^ergs^^, could be considered as 

reasonable values ( Willott et al. 1999{) . This may imply that if the beamed EC/CMB model for the 

255 large-scale jet in NGC 6251 with the parameters as given above is realistic at all, the outer parts of 

256 the outflow may provide some contribution to the observed 7-ray flux of 2FGL J 1629.4-1-8236, at 

257 least at the highest photon energies within the LAT range. 



253 Lp ~ 



258 



4.2. North- West Lobe 



259 Next we consider the case of the NW lobe of NGC 6251 likely associated with the Fermi- 

260 LAT source 2FGLJ1629. 4-1-8236. For the purpose of the spectral m odeling we measur ed radio 

261 fluxes for the lobe from published WSRT mdi^^ at 327 and 610 MHz jMack et al.lll997al) . Addi- 

262 tionally, we obtained a new 1.56 GHz map at 45" resolution with a total 5 hr integration 

263 from combining archival D-array observations obtained on 1985 December 1,5, and 1986 Jan 18 

264 (programs AB346 and TEST). Our measurements, S'327mhz = 3.10 Jy, ^eioMHz = 1-75 Jy, and 

265 5*1.56 GHz = 0.907 Jy, all assuming 10% errors, were measured using the same source extraction 

266 regions as in the Suzaku analysis (se e Figure 2) 0. The lobe also has a higher frequency radio flux 

267 measurement by Mack et al.' ( 1997a ). 5*10.55 ghz = 213.5 ± 26.4 mJy, and is detected and imaged 

268 down to 38 MHz (Rees,il990i) : due to the low resolution of the latter map, we summed the tabulated 



^Calibrated data sets obtained from the NRAO VLA Archive Survey. The National Radio Astronomy Observatory 
is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. 

"^We note that the radio source extraction region includes a relatively bright radio hotspot mentioned already in 
section 3.1. By means of careful analysis of the available radio maps we have estimated the contribution of this 
hotspot to the total radio emission of the lobe as 10 % at most. i.e. at the level of the assumed radio flux measurement 
uncertainty 
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269 fluxes for the 3 components listed west of the central radio source to derive an upper limit to the 

270 lobe flux, SssMRz < 39.6 Jy. Note that the emerging radio spectral index 0:0.326-10.55 ghz — 0.8 

271 is roughly consistent with the X-ray spectral index of the lobe measured by Suzaku . The lobe 

272 is approximated as a sphere with the radius = 185 kpc, so that the emission region volume is 

273 Ve = (4/3)7ri?f . 

274 For the lobe, we use a similar model as for the outer jet region discussed above, but include in 

275 addition extragalactic background light (EBL) photons within the infrared-to-optical range as tar- 

276 get photon field for the in verse-Compton scatter ing ('EC/(CMB+EBL)' model). In the calculations 



277 we use the EBL model of iMazin & Raud (|2007|) and account for Klein-Nishina effects. The choice 



278 of a particular EBL parameterization does not affect the model results because of the negligible 

279 contribution of the EC/EBL emission component to the 7-ray emission of the lobe within the LAT 

280 energy range. We neglect the contribution of the host galaxy to the target photon field because of 

281 the large distance of the analyzed NW lobe from the center (15' ~ 450 kpc). Finally, we assume 

282 homogeneous distributions of the radiating electrons and of the magnetic field within the lobe. 

283 In the framework of the applied simplified model, the 7-ray flux of 2FGL J1629. 4-1-8236 can 

284 be accounted for reasonably well (see Figure 4, thick gray line) by the inverse-Compton emission 

285 of the lobe for the following parameters: the equipartition ratio T]eq = Ue/Us = 45, magnetic field 

286 intensity B = 0.37 fiG, and the electron energy distribution of a double-broken power-law form 

287 rae(7) oc 7"''i for 7mm < 7 < 7br and Uei^y) oc 7""*' exp[-7/7max] for 7 > 7br with si = 2.0, 

288 S2 = 2.5, 7min = I7 7br = 3 X 10^, and 7jnax = 10^. With the above model parameters, the lobe 

289 pressure stored in ultrarelativistic electrons and magnetic field is Pc+B — 8.4 x 10^^^ dyncm^^. 

290 The total energy of the structure, being a sum of the work done in displacing a volume Vi of 

291 surrounding gas at pressure Pc+B, namely Pc+B Ve, and the energy of the material inside the cavity, 

292 Pc+B ye/il ~ 1) (assuming ultrarelativistic equation of state with the adiabatic index 7 = 4/3), is 

293 then i^tot — 4pc+B ~ 3 x 10^^ erg. The lobe's lifetime can be estimated as ti ~ Etot/Lj ~ 

294 1 X (Lj/10^^ erg s~^) Myr. We note that any significant proton contribution to the lobe pressure 

295 would increase the evaluated values of Etot and ti. 

296 The EC/(CMB-i-EBL) model fits the data quite well, and the emerging lobe's parameters seem 

297 reasonable. The implied departure from the energy equipartition is rather large, r^eq = 45, but still 

298 with in the range 1 — 100 establ i shed for the lobes in radio galaxies wit h X-ray measurements alone 



299 (e.g. iKataoka & Stawarzll2005l : ICroston et alj|2005l : llsobe et al.ll201 ih . Furthermore, the expected 



300 cooling timescale for the electrons emitting the observed > 1 GeV photons, i.e., electrons with 

301 energies 7 < 7max = 10^, is roughly tcooi — 2.1 x (7/10^) ^ Myr. Meanwhile, assuming a typical 

302 drift velocity of electrons in the lobe (or expansion velocity of the lobe) fiobc ^ 0.1 c, one can infer 

303 that these 7-ray emitting electrons can travel only d 65 x (t>iobc/0.1 c) x (7/10^) ^ kpc before 

304 they cool radiatively. This scale is much smaller than the extension of the NW lobe in NGC 625 1 
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305 (~ 0.5 Mpc). Interestingly, this finding resembles the case of the giant lobes in Centaurus A, and 

306 similarly suggests efficient in-situ re-accel eration of the radia ting particles within the whole vol- 

307 ume of the giant lobe (see the discussion in lAbdo et al.ll2010dh . Moreover, in both the Centaurus A 

308 and NGC 6251 radio galaxies, the 7-ray detections imply that the lobes emit one to two orders of 

309 magnitude more energy in 7-ray s than at radio/sub-mm wavelengths. 

310 Finally we note that the non-thermal pressure of the NW lobe in NGC 6251, ~ 10~^^ dyncm"^, 

311 seems comparable with/bit less than the thermal pressure of the surrounding gaseous medium, 

312 pth ~ (10~^^ — 1 0~^^) dyncrn"^ at t he distances 200" — 1000" from the core. The thermal pressure 

313 was estimated by Evans et al.l(|2005[) based on the extrapolation of the X-ray halo profile associated 

314 with the group of g alaxies including NGC 6251 and detected by Chandra and XMM-Newton (see 

315 also in this context iBirkinshaw & Worralllll993l : ISambruna et al.ll2004ah . Such an extrapolation 

316 may not provide a realistic estimate however, and should rather be considered as an upper limit 

317 for the ambient medium pressure. This, together with a likely contribution of (mildly) relativis- 

318 tic protons to the lobe's internal pressure, may suggest that either the extended radio structures 

319 in NGC 6251 are in a pressure equilibrium with the environment, or that these are over-pressured 

320 cavities, so that the expansion of the giant lobes in the system is still pressure-driven on Mpc scale 

321 (basically beyond the extension of the thermal halo of the surrounding group of galaxies). 



322 5. Conclusions 

323 We have presented a Suzaku X-ray observation of the NW lobe in a nearby radio galaxy 

324 NGC 6251. Through this observation, we found for the first time non-thermal diffuse X-ray emis- 

325 sion associated with the NW lobe. Since the error circle of the 7-ray source 2FGL J1629. 4-1-8236 

326 contains both the NW lobe and outer jet region, we discussed the possible origin of 7-ray emission 

327 by detailed modeling of the spectral energy distributions assuming a synchrotron plus inverse- 

328 Compton emission on CMB background photons. We argued that, at least at energies below 10 

329 GeV, the observed 7-ray emission is well explained by non-thermal emission from the NW lobe 

330 rather than from the outer jet region, with reasonable physical parameters (such as magnetic field, 

331 electron density etc) expected in radio lobes in general. Since the spatial extent of the NGC 6251 

332 radio lobe is quite large (~ 0.5 Mpc), and the expected radiative cooling time is much shorter 

333 for 7-ray emitting electrons, efficient in-situ acceleration of radiative particles is necessary within 

334 the whole volume of the giant lobe. We also briefly discussed the non-thermal pressure stored 

335 in ultrarelativistic electrons and magnetic field (not considering a likely contribution of relativistic 

336 proton) of the NW lobe in NGC 625 1 , which is estimated as ~ lO^^'^ dyn cm^^ from the model, and 

337 it seems comparable with or somewhat less than the thermal pressure of the surrounding gaseous 

338 medium. 
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Table 1: Siizcikii XIS observation log. 



Target Name 


R.A. [deg] 


Dec. [deg] 


Exposure [ks] 


Obs. Start Time (UT) 


NGC6251_LOBE 


246.5880 


82.6370 


34.2 


2011 Apr 15 06:37:27 


NGC6251_L0BE_BGD1 


247.8200 


82.9170 


18.8 


2011 Apr 16 00:38:47 


NGC625 1 iOBE_BGD2 


245.4460 


82.3530 


12.1 


2011 Apr 16 10:25:32 



Table 2: Fitting parameters for the power-law model. 



Parameter 


Value 


Nn [cm-2] 


5.54 X lO^o (fixed) 


r 


1.90 ±0.15 


i^ikev [keV-^ cm-2 s-1] 


(9.11 ±0.91) X 10"^ 


X^/d.o.f 


20.54/14 




1.14 X 10-^ 


-^0.4-7.5 kev [ergcm-2 s~^] 


4.07 X 10-^3 


-^0.4-7.5 kev [ergs~^] 


5.4 X 10^^ 
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Fig. 1. — Suzaku X-ray image of the NW lobe in NGC 6251 (vignetting and exposure corrections 
applied). Data from XIS and XIS 3 are summed in the 0.4 — lOkeV energy band. The image 
shows the relative excess of smoothed photon counts (arbitrary units indicated in the bottom bar) 
and displayed with square root scaling. The green contours denote th e large-scale radio structure 
of the source observed by WSRT at 55" resolution (|Mack et al.lll997al) and indicate levels of 8, 31, 
54, 77 and lOOmJy/beam. The red ellipse denotes the 95% position error of 2FGL J 1629.4-1-8236. 
The position of the radio core of NGC 6251 is marked by the black cross at the center and the 
adjacent 'outer jet region' is marked by the magenta ellipse. 
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Fig. 2. — Same as Figure 1, but zooming into the NW lobe region. The Suzaku source extraction 
region is denoted by the cyan contours. The point X-ray sources removed from the analysis are 
marked by cyan circles with red stripes. 
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Fig. 3. — Suzaku XIS spectra of the diffuse emission component coinciding with the NW lobe of 
NGC 6251 in the photon energy range 0.4 — 7.5 keV, fitted by the model wabs x power-law. 
The XISO + XIS3 and XISl spectra are shown in black and red, respectively. 
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Fig. 4. — Broad-band spectral energy distributions of the outer jet and the NW lobe in NGC 6251, 
including the Ferm/-LAT sou rce 2FGL J1629.4+832 6. The radio and X-ray data for the outer 
jet (blue bow-ti es) are frorn dSarnbruna et alJ l2004aL based on the maps presented in Perley et 
al. 1984), and Evans et al.1 (|2005h . respectively (see §4.1). The radio fluxes of the NW lobe 
(red squares and red arrow) are described in section 4.2. The X-ray fluxes for the lobe, as found 
in this paper, are given in the 0.5 - 1.0 keV, 1.0 - 2.0 keV, 2.0 - 4.0 keV, and 4.0 - 8.0 keV 
bins (open red circles). Finally, the LAT fluxes of 2FGL J 1629.4-1-8326 (magenta stars) are taken 
from the 2FGL catalog kbdo et al.ll2011f) . Thin solid line represents the beamed IC/CMB model 
for the outer jet (see §4.1), assuming the association with 2FGL J1629.4-I-8326. Thick gray line 
represent the IC/(CMB-i-EBL) model for the NW lobe (see § 4.2), again assuming the association 
with 2FGL J 1629.4-1-8326. 



